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FOREWORD 


This  manual  presents  a computer  program  used  at  the  Naval  Surface 
Weapons  Center,  Dahlgren,  for  the  purpose  of  studying  the  feasibility  of 
the  Wiener-Kolmogorof f filter  as  applied  to  radar  altimetry  data. 

The  task  was  performed  as  an  aid  to  a preliminary  study  in  simulating 
and  interpreting  radar  altimetry  data.  This  study  was  conducted  by  the 
Astronautics  Division  under  the  guidance  of  Hr.  Samuel  Smith  III,  Dr. 
Charles  Cohen,  and  Dr.  Bemd  Zondek. 

Financial  support  for  the  task  was  provided  by  an  ARPA  grant  for  Re- 
search and  Development  studies  of  satellite  radar  altimetry. 

Mr.  Robert  N.  Learn,  Mr.  Frank  V.  Werme,  Jr.,  and  Mr.  Samuel  L. 

Smith,  III  are  acknowledged  as  reviewers  of  this  manual. 
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S EC  I ION  1.  General  Description 

1.1  Purpose  of  the  Program  Maintenance  Manual.  The  objective  for  writing 
this  Program  Maintenance  Manual  for  "A  Satellite  Radar  Altimetry  Data 
Analysis  Computer  Program  for  Geodetic  Applications"  is  to  provide  the 
maintenance  programmer  personnel  with  the  information  necessary  to  effectively 
maintain  the  system. 

1.2  System  Application.  During  the  program  design  stage,  this  program 
was  called  the  "Radar  Data  Analysis  Program."  It  has  the  code  name  RADAP 
for  identification. 

Basically,  RADAP  is  used  to  develop  and  evaluate  analysis  techniques 
for  applying  the  Wiener-Kolmogorof f Filter  to  recover  deflections  of  the 
vertical  and  geoid  heights  from  radar  altimetry  data  at  sea. 

In  order  to  study  the  filter  feasibility,  RADAP  processes  synthetic 
data  generated  by  a computer  program,  SARAS,  which  is  described  in  Reference 
1.  This  controlled  synthetic  data,  when  processed  by  RADAP,  permits  one 
to  readily  evaluate  the  qualities  of  the  filter.  These  results  can  then 
be  used  to  optimize  the  filter  characteristics.  Reference  3 gives  an  account 
of  the  filter  evaluation  for  various  types  of  synthetic  altimetry  input 
data. 

1.3  Equipment  Environment.  The  computer  program  is  coded  for  the  CDC 
6700. 

1.4  Program  Environment.  Coding  of  the  computer  program  is  in  the  Ex- 
tended Fortran  IV  language  which  operates  with  the  Scope  3.3  system. 

1.5  Conventions. 

a.  The  code  name  for  this  program  is  RADAP. 

b.  This  program,  RADAP,  requires  an  input  file  as  generated  by  the 
computer  program  SARAS;  see  Reference  1.  However,  real  radar  altimetry 
data  can  also  be  processed  by  this  program  if  it  has  the  appropriate  input 
format. 


1 


SECTION  2. 


System  Description 


2.1  General  Description.  The  purpose  of  this  program  is  to  analyze  the 
simulated  radar  altimeter  data  produced  by  the  Satellite  Radar  Altimetry 
Simulation  (SARAS)  computer  program.  This  analysis  is  accomplished  by 
applying  a Wiener-Kolmogorof f filter  (developed  in  NWL  TR-2626;  "Accuracy 
of  Deflection  of  the  Vertical  Derived  from  Satellite  Altimetry"  by  Charles 
J.  Cohen  and  Bernd  Zondek,  for  the  TASC  autocorrelation  function)  to  the 
noisy  simulated  data  and  evaluating  the  accuracy  of  the  geoid  heights  and 
the  deflections  of  the  vertical  obtained  against  the  "true"  values.  If 
this  program  performs  well,  it  will  form  the  heart  of  an  operational  system 
to  analyze  data  produced  by  the  GEOS-C  and  other  satellites. 

2.1.1  Program  Description.  Mathematical  Formulation 


2. 1.1.1  Input. 

SARAS  Output  Tape  (File) 

x(i)  - Data  Point  Position  (km) 

Y(i)  - True  Geoid  Height  (m) 

YN(i)  - Noisy  Geoid  Height  (m) 

D(i)  - True  Deflection  of  Vertical  (arc  sec) 

N - Number  of  Data  Points 

2. 1.1.2  Controls. 

a.  Select  Geoid  Height  and/or  Deflection  of  Vertical 

b.  Filter  Fit  Span  (±L) 

(number  of  points) 

2. 1.1.3  Constants. 

0 - Reciprocal  decay  distance  (km-1) 

- Standard  deviation  of  single  point  distance  measurement  (meters) 

0n,  - Standard  deviation  of  deflection  of  the  vertical  (arc  sec) 

^s  - Filter  fit  spacing  (km)  - must  be  an  integral  multiple  (£)  of 
data  spacing 

1 - as/ax  (integer) 


2. 1.1.4  Precomputations. 


f 16-ns 
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(1) 
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(2) 

°1 

a/2i  [2  jr^ 

R + R2  + (R  - 

2)]‘ 

(3) 

Ps  * 

a/3  [2  yr^ 

R 4 

- R2  - (R  - 
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j 
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2. 1.1. 5 Geoid  Heights  Computation. 

Compute  Filter  Weights , (x(  ) 
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Compute  Filter  Weights,  WN»  (x^ ) 


(9)  F'(rj)=  | ^(Cos  Plr3  - /Tsin  Pir3)  e 


•OnT 


1+R  r 


- e 


] 


(10)  W^'CSj)  = sgn(s^)  g2F  (r ^ ) 
Plot:  / ( s 5 ) vs  Xj  , -L^s  ^ Sj  <;  L^s 


Title  for  Plots:  Weighting  Function  for  Deflection  of  Vertical 

Run  No.  , Date  


Data/Noise  Type 


-/■ 


Compute  Smoothed  Deflection  of  Vertical,  N (x. ) 

For  each  x{  from  i = -(N  - l)/2  + (L  to  i = (N  ■ l)/2  — 

A ✓ 

compute  the  estimated  deflection  of  the  vertical,  N (x^),  by: 


a / 


(11)  N (xt)  = AS  r WN'(Sj)  YN(i  + JA) 

J = 'L 


2. 1.1. 7 Statistical  Comparisons. 

A 

If  Geoid  Heights:  s(  = N(x(  ) 


s'  = V(i) 


If  Deflection  of  the  Vertical: 


Sj  = N (x,  ) 


s,'  - D(i) 


1 


i 


Compute  residuals  oC  "Estimated"  and  "True"  values  by 

(12)  Ast  = Sj  - s/ 

Print  the  Maximum  Value  of  Residual,  [Max  As;  j 
Compute  Average  Value  of  Residuals,  As,  by 


(13)  As  = 1/  (N  - 2/L)  ASj 

5 * 1 


Print  the  Average  Value  of  Residual,  As 
Compute  Standard  Deviation  of  Residuals,  crs,  by 


(14)  - (I/O.  - 2lL)  (4S_  . -)s) 


Print  Standard  Deviation  of  Residuals, 


5 


2.1.2  Program  Description.  Definitions  and  Structure 

2.1.2. 1 General  Structure.  The  computer  program  RADAP  consists  of  an 
executive  routine  called  RADAP  and  seven  subroutines.  These  subroutines 
are  called  TLIST,  GHCMPTS , CMPSGH,  STATC0M,  DVCMPS,  PL0T,  and  SCALE.  A 
chart  showing  the  interrelationships  between  the  routines  is  given  in 
F i gure  1 . 


INPUT 

CARDS 


RADAP 


(FILE 

V (2) 


GHCMPTS 

(2) 


CMPSGH 

(3) 


STATC0M 

(4) 


DVCMPS 

(5) 


CMPSGH 

(6) 


/ PRINT 

/g.h.w. 


I PRINT 
G.H.E. 


f PRINT 

G.H.S. 


! PRINT 
V.D.W. 


I PRINT 
V.D.E. 


/ PRINT 

V.D.S. 


De  finitions: 


File  1 

File  2 

RADAP 

TLIST 

GHCMPTS 

CMPSGH 

STATCOM 

DVCMPS 

PLOT 

SCALE 


- Input  file  to  RADAP  containing  a header  record  and  two  data 
records.  The  first  data  record  contains  the  along  track  dis- 
tances and  the  noisy  geoid  heights.  The  second  data  record 
contains  the  smooth  geoid  heights  and  the  true  vertical  de- 
flections. 

- The  output  file  from  RADAP  containing  four  records.  The  first 
two  are  header  records  and  the  last  two  records  are  data  re- 
cords. Records  three  and  four  contain  the  geoid  height  and  the 
vertical  deflection  residuals,  respectively. 

- Executive  or  control  routine  for  the  program. 

- Reads  the  data  records  of  File  1 into  the  appropriate  array 
locations.  These  data  records  contain  the  Along  Track  Distance 
Array,  the  Noisy  Geoid  Height  Array,  the  Smooth  Geoid  Heights 
Array,  and  the  Vertical  Deflection  Array. 

- Computes  the  geoid  height  filter  weights,  (s ^ ) , for  all  data 
points  and  plots  these  weights  vs  along  track  distance  xj . 

- Computes  the  smoothed  Geoid  Height,  N^),  or  the  smoothed  De- 
flection of  the  Vertical,  ft'Cx^,  arrays  dependent  upon  which 
set  of  filter  weights,  WN(x  ) or  W^CXj),  are  residing  in  the 
filter  weight  array  (WN). 

- Residual  statistics  for  Geoid  Height  or  Deflection  of  the 
Vertical  are  computed  in  terms  of  maximum  average  and  standard 
deviation  values. 

- Computes  the  vertical  deflection  filter  weights  WN/(s  ),  for 
all  data  points  and  plots  these  weights  vs  along  traclc  dis- 
tances xf  . 

- A control  routine  for  all  RADAP  plotting  requirements.  It  con- 
trols the  titling,  labeling  of  axes,  and  plotting  of  points  for 
each  graph. 

- Scales  the  ordinate  values  of  each  plot  so  that  the  maximum 
ordinate  has  the  value  0.5  or  1.0. 
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According  to  Figure  1.,  the  executive  routine  RADAP  reads  input  cards, 
reads  file  L,  writes  on  file  2,  and  calls  the  seven  subroutines.  These 
seven  subroutines  are  called  in  the  order  as  numbered  in  Figure  1. 

Subroutine  CMPSGH  is  called  twice,  once  to  compute  and  print  the 
geoid  height  estimates  (G.H.E.)  and  once  to  process  the  vertical  deflection 
estimates  (V.D.E.). 

Subroutine  STATC0M  is  also  called  twice,  once  to  process  geoid  height 
residuals  and  once  to  process  vertical  deflection  residuals.  The  printout 
by  STATC0M  consists  of  the  residual  statistics  maximum,  average,  and  standard 
deviation  (G.H.S.  and  V.D.S.). 

Printing  is  also  done  by  subroutines  GHCMPTS,  DVCMPS,  and  STATC^M. 
GHCMPTS  prints  the  geoid  height  weights  (G.H.W.).  DVCMPS  prints  out  the 
vertical  deflection  weights  (V.D.W.). 

A plotting  subroutine,  PL0T,  is  called  by  subroutines  GHCMPTS,  STATC0M, 
and  DVCMPS.  Subroutine  SCALE  is  called  by  the  PL0T  subroutine  in  order  to 
establish  a suitable  graphical  scale.  Four  graphical  plots  are  produced 
by  the  program. 

2.2  Detailed  Description.  Program  RADAP  will  now  be  described  in  detail. 

To  do  this,  a narrative  description  will  be  given  for  the  program  routines 
RADAP,  TLIST,  CHCMPTS , CMPSGH,  STATC0M,  DVCMPS,  PL0T,  and  SCALE. 

Within  the  narrative  description,  use  is  made  of  the  routine  variables 
as  given  in  Appendix  A.  The  narratives  follow  the  macro-flowcharts  for 
each  routine  as  given  in  Appendix  B. 

Further  information  for  each  program  routine  is  contained  in  the  Ex- 
tended Fortran  IV  program  listings  of  Appendix  C.  Appendix  D gives  a 
sample  set  of  program  output.  The  output  will  be  discussed  more  fully  in 
Section  3. 

2.2.1  RADAP  Narrative  Description. 

a.  This  routine  is  the  executive  or  control  routine  for  the  program 
flow  during  execution. 

b.  The  routine  functions  as  follows:  At  the  start,  the  routine  prints 

out  three  lines  for  program  identification  purposes.  This  identification 
defines  the  program  by  name,  version  number,  and  a title  indicator  for  card 
data  input.  Next,  the  routine  reads  and  prints  the  contents  of  the  two 
input  cards  which  contain  radar  analysis  constants  and  controls. 

Initialization  of  routine  constants  takes  place  and  file  (1)  is  re- 
wound and  read  for  the  function  selector  value  NFS.  NFS  can  take  on  the 
values  1 through  5 and  represents  a geological  feature  simulated  in  the  SARAS 
program.  The  NFS  values  are  identified  as  follows: 
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NFS  = 1 
NFS  = 2 
NFS  = 3 
NFS  = 4 
NFS  = 5 


symmetric  ocean  trough  or 

escarpment 

inclined  plane 

sea  waves,  autocorrelated 

sea  waves,  autocorrelated 


crest 


- Zondek  method 

- Cohen  method 


The  geological  feature  names  require  2 or  3 words  of  Hollerith  code 
to  represent  these  titles  and  the  printout  of  these  names  requires  dif- 
ferent formats.  Files  (1)  and  (2)  are  rewound.  A test  is  now  performed 
on  the  NFS  value.  If  NFS  is  a 1 or  2,  then  IT  = 3;  where  IT  is  the  num- 
ber of  words  describing  the  geological  feature  being  processed.  A test 
on  IT  will  cause  the  program  to  go  to  the  appropriate  set  of  I/O  statements 
to  read  and  print  the  header  record  from  file  (1)  that  is  associated  with 
the  geological  feature  defined  by  the  value  of  NFS. 


A test  of  the  index  IC,  the  option  control  for  producing  an  output 
file,  is  made  after  the  header  record  of  file  (1)  is  read.  If  IC  = 0, 
then  no  output  file  is  produced.  Otherwise,  the  header  record  of  file  (1) 
is  written  on  file  (2)  and  it  is  also  printed. 


II,  the  total  number  of  data  points  minus  the  number  of  filter  fit 
points,  is  computed  and  is  written  on  file  (2)  along  with  other  values  ob- 
tained from  input  card  data.  This  written  record  on  file  (2)  is  the 
second  header  record  and  this  information  is  also  printed. 


At  statement  250,  the  RADAP  routine  calls  the  program  subroutines  in 
the  order  which  follows: 


TLIST , GHCMPTS , CMPSGH , 

STATC0M, DVCMPS , CMPSGH , 

STATC0M. 

After  all  appropriate  subroutines  have  been  called  and  file  (2)  has 
been  rewound,  the  library  subroutines  CRTID  and  EXITG  are  called  to 
identify  the  graphical  output  and  close  out  the  plotting  procedures.  The 
program  is  then  terminated. 

c.  The  storage  required  by  this  routine  is  distributed  as  follows: 

Program  length  742g  (482  ) 

Buffer  length  12132ft  (5210  ) 

Common  length  112202fi  ( 38018  ^ 0 ) 


2.2.2  TLIST  Narrative  Description. 

a.  TLIST  reads  the  input  data  records  2 and  3 from  file  (1)  into 
appropriate  arrays  for  further  program  processing. 

1.  Record  2 consists  of  the  x and  YN  arrays,  where  x = Data 
Point  Positions  (km)  and  YN  = Nosiy  Geoid  Heights  (m) . 
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2.  Record  3 consists  of  the  Y and  D arrays  where  Y = True  Geoid 
Heights  (m)  and  D .=  True  Deflections  of  Vertical  (arc  sec). 

b.  After  records  2 and  3 are  read  into  the  appropriate  arrays,  file 
(1)  is  rewound  and  control  returns  to  the  executive  routine. 

c.  This  subroutine  has  a storage  length  of  51(_  or  41iq  words. 

2.2.3  GHCMPTS  Narrative  Description. 

a.  The  purpose  of  this  subroutine  is  to  compute  and  control  the 
plotting  of  the  geoid  height  filter  weights. 

b.  Initially,  the  subroutine  calculates  and  sets  local  program  con- 
stants. The  standard  deviation  of  deflection  (°n0  is  converted  from  arc 
sec  to  radian  measure.  Standard  deviation  (CT3)  for  a single  point  distance 
in  meters  is  converted  to  kilometers.  Next,  the  geoid  height  filter 
weights  are  computed  and  stored  in  array  WN  by  a loop  from  1 to  LL.  LL 

is  equal  to  the  number  of  filter  fit  points  plus  one. 

These  filter  weight  values  are  printed  as  part  of  the  hardcopy 
output.  Next,  a call  is  made  to  the  PL0T  subroutine  for  plotting  of  the 
geoid  height  filter  weights  versus  the  along  track  distance.  The  program 
processing  then  returns  to  RADAP. 

e.  The  program  length  is  244  (164  ) and  common  storage  is  111670 

(37816JO). 

2.2.4  CMPSGH  Narrative  Description. 

a.  This  subroutine  computes  the  geoid  height  or  vertical  deflection 
estimates  N^)  or  N(x;  ) depending  on  which  one  of  the  weights  are  in  the 
WN  array. 

b.  The  argument  list  to  this  subroutine  are  defined  as  follows: 

WN  = array  storage  for  filter  weights 
YNT  = array  storage  for  calculated  geoid  heights  or  vertical 
deflections 

NP  = two  times  number  of  filter  fit  points  plus  one 
YN  - array  storage  for  noisy  geoid  heights 
IC  = 1,  CMPSGH  prints  geoid  height  estimates; 

= 2,  CMPSGH  prints  vertical  deflection  estimates. 

The  first  part  of  this  subroutine  initializes  and  computes  local 
and  global  program  constants. 

Next,  the  counter  J and  the  summation  variable  YNN  are  set.  Com- 
putation is  accomplished  by  a nested  loop  procedure.  The  inner  loop 
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evaluates  the  appropriate  sum  for  YNN.  A continuation  in  the  outer  loop 
computes  the  geoid  height  or  vertical  deflection  estimates  and  stores  them 
in  array  YNT. 

Attest  on  the  index  IC  is  made  and  the  appropriate  estimates 
N(x  ) or  N(xj)  are  printed  and  the  subroutine  ends.  The  printout  is  done 
with  the  CDC  computer  system  library  routine  PDUMP.  Each  estimated  value 
has  an  associated  storage  location  value  printed  also. 

c.  Subroutine  storage  is  142e(98l0)  and  common  length  is  20p(16iQ). 

2.2.5  STATC0M  Narrative  Description. 

a.  This  subroutine  computes  the  residual  statistics  for  geoid  heights 

or  vertical  deflections.  The  statistics  computed  for  a set  are  *-he 

residual  maximum,  average,  and  standard  deviation;  designated  respectively 
by  | MAX  asi  ] » A®>  and  r • 

b.  First,  a test  is  made  on  the  index  IC.  If  IC  = 1,  the  geoid 
height  residuals  are  computed.  IC  = 2 means  that  the  vertical  deflection 
residuals  are  computed.  A test  of  index  IIC  is  made  if  the  residuals  are 
to  be  written  on  the  output  file  (2).  IIC  ^ 0 implies  that  the  output 
file  is  written.  This  test  is  made  for  both  types  of  residuals. 

The  next  operation  involves  finding  the  maximum  absolute  value  of 
the  residuals.  A calculation  for  the  average  value  of  the  residuals  is 
then  performed  followed  by  a calculation  of  the  residual  standard  deviation. 

Another  test  of  index  IC  is  made  for  printing  and  plotting  the 
appropriate  residuals  and  residual  statistics.  STATC0M  processes  only  one 
type  of  residuals  each  time  it  is  called  from  RADAP. 

c.  Subroutine  length  is  301  (193  ) and  common  length  is  111670 

(37816l0). 

2.2.6  DVCMPS  Narrative  Description. 

a.  The  purpose  of  this  subroutine  is  to  compute  and  control  the 
plotting  of  the  vertical  deflection  filter  weights. 

b.  DVCMPS  has  the  same  sequence  of  events  as  subroutine  C.HCMPTS. 

See  the  narrative  description  of  GHCMPTS. 

c.  Subroutine  storage  required  is  301  (19313)  for  program  length 
and  111670r (37816]  ) for  common  length. 


2.2.7  PLOT  Narrative  Description. 

a.  This  subroutine  acts  as  the  control  for  writing  the  title, 
labeling  the  axes,  and  plotting  the  points  for  each  graph.  Details  of 
the  graphical  subroutines  used  by  PLOT  are  contained  in  Reference  4.  The 
subroutine  argument  list  is  defined  as  follows: 

ZN  = ordinate  array  for  storage  of  true  geoid  heights,  true 

vertical  deflections,  weighted  geoid  heights,  or  weighted 
vertical  deflections 
X = data  point  position  abscissa  array 

NL  = number  of  data  points  to  be  plotted 

LT  = graph  number 

LT1  = number  of  plots  per  graph  (1  or  2). 

NFS  = selector  for  geological  feature 

LK  = location  for  GRF  argument  M. 

The  PLOT  subroutine  has  the  capability  of  plotting  any  of  the 
following  ordinate  value  arrays  versus  along  track  distance  (x; ) : 

(a)  weighted  geoid  height 

(b)  weighted  vertical  deflection 

(c)  true  geoid  height  and  geoid  height  residuals  (on  the  same 
graph) 

(d)  true  vertical  deflection  and  vertical  deflection  residuals 
(on  the  same  graph) 

However,  only  one  of  the  above  are  processed  at  each  call  to 
PL0T.  Titling  of  plots  is  dependent  upon  the  following  values  of  plot 
number  LT: 

1 - weighting  function  for  geoid  height 

2 - weighting  function  for  vertical  deflection 

3 - true  geoid  height  and  geoid  height  residuals 

4 - true  vertical  deflection  and  vertical  deflection  residuals 

b.  The  first  executed  statement  in  the  subroutine  is  a call  to 
the  library  subroutine  MODESG.  This  initial  call  to  MODESG  sets  the  Z 
array  to  the  default  values  as  used  by  the  plotting  subroutines. 

A test  is  performed  on  the  value  of  LT.  LT  > 1 implies  that  this 
is  the  second  pass  through  the  subroutine  PLOT  and  the  plot  identification 
procedures  are  by-passed.  If  LT  = 0 or  1,  then  it  is  the  first  pass  through 
PLOT  and  the  film  identification  is  read  from  an  input  card.  A call  is 
now  made  to  the  film  identification  subroutine  CRTID.  Another  input  card 
is  read  which  contains  thirty  alphanumeric  characters  which  define  the 
graph  title,  noise  type  and  standard  deviation,  and  the  plot  interval. 

Next,  the  maximum  absolute  ordinate  value  is  determined  from 
ordinate  array  ZN.  Index  LT1  is  now  tested.  When  LT1  = 1,  then  there  is 
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only  one  plot  per  graph.  If  LT1  = 2,  then  there  are  two  plots  per  graph 
and  the  maximum  absolute  ordinate  value  must  be  determined  from  both  sets 
of  data;  the  true  values  and  their  residuals  (ZN  and  WN  arrays). 

The  maximum  and  minimum  abscissa  values  are  now  found.  And  next, 
the  SCALE  subroutine  is  called  so  that  the  scale  factor  KE  can  be  determined 
so  that  a new  convenient  maximum  ordinate  of  0.5  x 10*  or  0.1  x 10* , where 
x is  the  appropriate  exponential  power,  is  defined. 

The  graphical  axis  origin  is  now  set  in  locations  X0R  and  Y0R. 

DX  and  DY  are  set  for  the  abscissa  and  ordinate  tick  mark  spacings.  Scale 
factor  KE  is  now  stored  in  a floating  point  location  R. 

Another  test  is  now  performed  on  the  index  LT1  for  scaling  pur- 
poses. If  LTl  ■ 2,  then  there  are  2 plots/graph.  The  residuals  and  the 
ordinate  values  contained  in  the  ZN  array  are  scaled  and  stored  in  the 
WN  and  ZW  arrays,  respectively.  If  LTl  = 1,  only  the  ZN  array  is  scaled 
for  plotting.  Other  plot  control  values  are  also  initialized  before  setting 
up  the  graphical  subject  and  object  spaces. 

At  statement  300,  the  subject  space  is  established  in  raster  units 
and  the  object  space  is  established  in  normalized  units.  For  plot  control, 
calls  to  SETSMG  are  necessary  for  updating  the  Z array  before  calls  to  the 
subroutine  TITLEG.  Calls  to  NUMBRG  and  LEGNDG  are  made  in  order  to  print 
the  scale  factor  on  the  graph  at  a particular  location.  The  control  vari- 
able LT  selects  the  appropriate  graphical  plot  labeling  depending  on  the 
number  of  the  graph  being  plotted. 

After  plotting  one  of  the  graphs,  the  processing  goes  to  statement 
700.  Here,  the  object  space  is  reset  and  a call  is  made  to  GRF.  Library 
subroutine  GRF  puts  on  the  axis  and  plots  one  curve  on  the  graph. 

If  LTl  = 1,  then  the  processing  skips  a call  to  GRF  since  no 
residuals  are  plotted.  After  statement  750,  a call  to  library  subroutine 
PAGEG  (skips  to  next  graphical  page)  is  made  and  the  subroutine  PL0T  returns 
control  to  RADAP. 

c.  Storage  needs  are  16015^  (718110)  words  for  program  length  and 
112200^  (3801610)  words  for  common  length. 

2.2.8  SCALE  Narrative  Description. 

a.  This  subroutine  is  used  to  scale  the  ordinate  values  to  be  plotted 
so  that  the  maximum  ordinate  has  the  value  0.5  or  1.0  appearing  on  the 
graph.  The  appropriate  exponent  is  also  determined  for  labeling  of  the 
graphical  ordinate  scale. 
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b.  Initially  the  variable  I,  used  to  store  the  integer  exponent  of 
the  maximum  ordinate  YMA,  is  set  to  zero.  The  decimal  part  of  YMA  is  set 
between  -1.0  and  0.0. 

The  exponential  part  of  YMA  is  stored  into  location  I.  If  the 
decimal  part  of  YMA  is  between  -0.5  and  0.0,  then  the  graphical  ordinate 
value  YMI  is  set  equal  to  -0.5.  Otherwise,  YMI  = -1.0  and  the  program 
processing  returns  to  subroutine  PL0T. 

c.  Storage  subroutine  length  is  63s(5110)  words. 

2.3  Processing  Time. 

The  computer  time  necessary  to  process  a total  of  6251  data  points 
is  164  central  processor  seconds  (CPS).  This  time  estimate  includes  the 
time  for  the  generation  of  an  output  file. 

Compilation  time  for  this  program  is  14.0  CPS. 


SECTION  3.  INPUT/OUTPUT  DESCRIPTION 


3.1  General  Description.  This  program  utilizes  input  data  from  a 
permanent  file  and  IBM  cards.  The  permanent  file  is  generated  by 
the  SARAS  program  described  in  Reference  1.  In  addition,  the  RADAP 
program  can  now  process  real  satellite  radar  altimetry  data  if  it  is 
in  the  appropriate  input  format. 

Input  card  data  contains  information  about  simulated  satellite 
radar  altimeter  observations.  Card  data  also  contains  design  values 
for  the  Weiner-Kolmogorof f filter. 

RADAP  output  consists  of  hard  copy  printout,  four  graphs,  and 
an  output  file,  on  option. 

3.2  Characteristics,  Organization,  and  Detailed  Description  of 
System  Data. 

3.2.1  General  Characteristics.  During  program  execution,  the  input 
to  the  program  consists  of  80  column  IBM  cards  and  a permanent  file. 
Output  of  the  program  includes  hard  copy  printout,  a permanent  file 
(on  option),  and  four  separate  sheets  of  graphical  output. 

a.  Routine  RADAP  first  prints  out  the  program  identification. 

It  then  reads  and  prints  out  the  input  card  data.  It  also  reads  and 
prints  out  the  header  record  from  the  input  permanent  file,  SARAS  data. 
If  the  output  permanent  file  option  is  on,  then  RADAP  writes  and  prints 
out  the  first  two  header  records. 

Subroutine  TLIST  reads  the  remaining  two  data  records  from 
the  input  data  file  into  appropriate  arrays  for  further  processing. 

Subroutine  GHCMPTS  prints  and  sets  up  the  call  to  plot  the 
geoid  height  weights. 

CMPSGH  prints  the  hard  copy  of  geoid  height  estimates  or 
vertical  deflection  estimates. 

STATC0M  prints  out  the  statistics  for  the  residuals  and 
initiates  a call  to  plot  true  values  and  residuals  for  either  geoid 
heights  or  vertical  deflections. 

DVCMPS  prints  and  initiates  the  call  to  plot  vertical 
deflection  weights. 

Subroutine  PL0T  reads  the  last  two  input  data  cards  and 
actually  calls  the  graphical  plotting  library  subroutines. 
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b.  Input  data  cards  convey  information  to  the  program  RADAP 
of  the  following  types: 

1.  Program  constants  needed  to  evaluate  the  filter. 

2.  Control  option  for  writing  permanent  output  file. 

3.  Film  identification,  data  for  graph  title,  and  a plot 
interval  numerical  value. 

The  input  permanent  file,  SARAS  data,  contains  simulated  satellite 
altimeter  observations  to  be  processed  by  the  program.  Hard  copy 
printout  and  graphical  output  is  used  for  both  monitoring  and 
evaluating  the  filter  qualities. 

A permanent  output  file  which  contains  the  geoid  height  and 
vertical  deflection  residuals  is  used  by  other  computer  programs  for 
further  analysis  (e.g.  harmonic  analysis). 

c.  Input  data  to  the  program  is  of  a fixed  or  static  type. 

The  program  does  not  alter  or  modify  this  data  in  any  way. 

d.  Storage  requirements  for  the  input  and  output  files  ate  about 
38000g  words. 

e.  The  input  permanent  file  data  is  limited  to  43000g  words 
since  the  program  loading  and  execution  uses  211000g  words. 

3.2.2  Organization  and  Detailed  Description.  The  internal  structure 

for  each  input/output  component  will  now  be  discussed.  In  addition,  l 

a set  of  sample  output  given  in  Appendix  D will  be  described. 

Input  data  cards  are  described  as  follows: 

Card  1:  Program  constants  read  by  the  RADAP  routine. 


Column 

Variable 

Description 

Format 

1-25 

BRDC 

Inverse  Correlation 
Distance  (1/km) 

E25.14 

26-50 

SDD 

Single  Pulse  Radar 
Standard  Deviation  (m) 

E25.14 

51-55 

NFS 

Function  Selector 

15 

56-60 

LSX 

Filter  Fit/Data  Spacing 

15 

61-65 

IC 

Output  File  Option 
0=no  file,  l=output  file 

15 
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Card  2:  Program  constants  read  by  the  RADAP  routine. 
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Column 

Variable 

Description 

Format 

1-25 

SDV 

Deflection  Standard 
Deviation  (arcsec) 

E25.14 

26-50 

FDS 

Data  Spacing  (km) 

E25.14 

51-60 

LFSP 

\ Number  of  Filter 
Fit  Points 

110 

Card  3: 

Film  identification  read  by  subroutine  PL0T 

Column 

Variable 

Description 

Format 

1-20 

H 

Film  Identification: 

NWXY  - columns  1-4 
XY  is  programmer  number 
XXXX  - columns  5-8 
program  name 
XXXX  - columns  9-12 
office  room  number 
XXXXXXXX  - columns  13-20 
calendar  date 

2A10 

Card  4: 

Graph  title 
subroutine 

data  and  plot  interval  as  read 
PL0T 

by  jJ 

Column 

Variable 

Description 

Format 

1-30 

DA 

Thirty  Alphanumeric 
Characters  for  Graph  Title 

3A10 

31-60 

TN 

Noise  Type  and  Standard 
Deviation  (m) 

3A10  ' 

61-65 

IP 

Plot  Interval,  i.e. 
Plots  Every  IP^h  Point 
From  the  Data  Set 

|S 

1 j 

The  structure  for  the  input  permanent  file,  SARAS  data,  can  be 
found  in  Reference  1,  page  B-8. 

An  output  file,  in  binary  form,  is  produced  if  the  output  control 
option  IC-1.  The  internal  structure  of  the  output  file  is  as  follows: 


if 


' 
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Record  1.  Header  Record  (1) 


Word  # 

Word  Content 

Word 

Format 

Description 

1 

NFS 

I 

Function  selector 

2 

NR 

I 

Run  number 

3 

DATE 

A 

Calendar  date 

4-5  or 

B(I) 

A 

Geodetic  feature  (number  of 

4-6 

words  depends  on  NFS) 

6 or  7 

NGS 

I 

Noise  type  generator 

7 or  8 

SN 

F 

Noise  standard  deviation 

8 or  9 

N 

I 

Total  number  of  data  points 

Record  2.  Header  Record  (2) 

Word  # 

Word  Content 

Word 

Format 

Description 

1 

NFS 

I 

Function  selector 

2 

BRDC 

F 

Inverse  correlation  distance 

3 

SDD 

F 

Single  pulse  radar  standard 
deviation 

4 

LSX 

I 

Filter  fit/Data  spacing 

5 

SDV 

F 

Deflection  standard  deviation 

6 

FDS 

F 

Data  spacing 

7 

LFS 

I 

Number  of  filter  fit  points 

8 

II 

I 

Total  number  of  data  points 
minus  total  number  of  filter 
fit  points 

Record  3.  Data  Record 

Word  # 

Word  Content 

Word 

Format 

Description 

1 

NP 

I 

Number  of  data  points 

2 

WN 

F 

Geoid  height  residual  array 

3 

WN 

F 

NP  + 1 

WN 

F 
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Record  4. 


Data  Record 


This  record  has  the  same  structure  as  record  3 except  that  the  array- 
values  are  vertical  deflection  residuals. 

The  hard  copy  print  out  and  the  graphical  output  of  the  program  will 
now  be  described.  A sample  of  program  output  is  contained  in  Appendix  D. 

On  page  1 of  Appendix  D appears  the  program  identification.  This  is 
printed  from  the  RADAP  routine  and  includes  the  program  name,  version,  and 
date.  At  the  top  of  the  next  page  appears  the  title  of  the  output  data. 

The  next  9 lines  gives  the  input  card  data  as  printed  from  the  RADAP 
routine. 

A printout  of  the  header  record  as  read  from  the  input  file  is  con- 
tained in  the  next  two  lines.  Following  in  the  next  3 lines  is  the  header 
records  for  the  output  file.  Both  input  and  output  file  header  records 
are  printed  from  routine  RADAP. 

The  next  set  of  data  is  printed  from  Subroutine  GHCMPTS.  This  print- 
out contains  the  values  for  geoid  height  weights  and  are  read  across  the 
page  for  each  line.  This  is  a PDUMP  which  begins  at  octal  location 
061260e  and  ends  at  0632408. 

The  next  2 lines  of  printout  are  input  card  data  as  printed  from  Sub- 
routine PL0T.  Each  line  of  values  are  identifiable  from  the  printout. 

A PDUMP  for  the  geoid  height  estimates  are  printed  from  Subroutine 
CMPSGH  beginning  at  octal  location  0755148  and  ending  at  103000e . 

Geoid  height  residual  statistics  appear  on  the  next  3 lines  and  are 
clearly  identified  by  titles.  These  values  are  printed  from  Subroutine 
STATC0M. 

Now,  the  same  PDUMP  printout  follows  for  the  vertical  deflections  as 
given  above  for  the  geoid  heights.  Subroutine  DVCMPS  prints  out  the 
vertical  deflection  weights.  The  vertical  deflection  estimates  are  printed 
from  Subroutine  CMPSGH.  STATC0M  now  prints  out  the  vertical  deflection 
residual  statistics. 

The  last  two  lines  of  printout  are  provided  by  the  system  library 
graphical  routines.  This  information  informs  the  programmer  as  to  the 
status  of  the  graphical  output.  A detail  explanation  of  this  printout 
can  be  found  in  Reference  4,  pages  4-50  to  4-53. 
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After  the  printout,  the  graphical  output  appears  as  shown  in  Appendix 
D.  Four  graphical  plots  are  produced.  The  first  two  graphs  are  plots  of 
geoid  height  and  vertical  deflection  weights  versus  the  along  track  distance. 

Graph  number  three  contains  two  separate  plots.  One  curve  is  a plot 
of  true  geoid  heights  versus  along  track  distance.  The  other  curve  on 
the  same  ordinate  scale  is  a plot  of  the  geoid  height  residuals  versus 
along  track  distance. 

The  fourth  graph  is  similar  to  the  third  graph  except  that  the  ordinate 
values  apply  to  true  deflection  of  vertical  and  deflection  of  the  vertical 
residuals . 
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Appendix  A 

Definition  of  Program  Variables 
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RADAP 


B - An  array  which  contains  the  name  of  the  geological  feature  being  pro- 
cessed. 

BRDC  - Inverse  correlation  distance  (1/km) 

D - An  array  used  to  store  the  true  vertical  deflections 
DATE  - Calendar  date  of  input  file  data 
FDS  - Data  spacing  (km) 

I - Do  loop  integer  variable 

IC  - Output  file  option  control  index 
IGF  - Lower  do  loop  index  variable 

II  - Numerical  value  of  total  number  of  data  points  minus  number  of  filter 

fit  points 

IT  - Upper  do  loop  index  variable 
L = LFSP  - k number  of  filter  fit  points 
LFS  - Number  of  points  for  filter  fit  plus  1 
LSX  - One  filter  fit/data  spacing 
N - Total  number  of  data  points 

NFS  - Variable  which  depends  on  the  geological  feature  being  processed 
(see  reference  1,  page  B-8) 

NGS  - Noise  type  generator  (see  reference  1,  page  B-8) 

NR  - Input  file  run  number  (see  reference  1,  page  B-8) 

SDD  - Single  pulse  radar  standard  deviation  (meters) 

SDV  - Deflection  standard  deviation  (arc  sec) 

SN  - Noise  standard  deviation  (see  reference  1,  page  B-8) 
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D - Array  for  smooth  vertical  deflection  values 
X - Array  for  along  track  distance  values 
Y - Array  for  smooth  geoid  height  values 
YN  - Array  for  noisy  geoid  height  values 
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AR  - Temporary  storage  variable 
BRDCQ  - Cube  of  inverse  correlation  distance 
BRDCS  - Square  of  inverse  correlation  distance 
DR  - Temporary  storage  variable 

F - Temporary  storage  variable  used  to  calculate  filter  weight 

FDS1  - Temporary  storage  of  a count  value  (J)  times  data  spacing; 

-LFSP  s J < LFSP. 

I - Do  loop  variable 

ID  - Variable  used  to  store  the  number  of  data  points  taken  up  by  number 
of  filter  points 

IL  - Total  number  of  data  points  on  one  side  of  middle  minus  number  of 
points  taken  up  by  filter  fit  points  plus  1. 

J - Variable  counter,  -LFSP  ^ J ^ LFSP. 

LL  - Two  times  number  of  filter  fit  points,  plus  one 

LM  - Number  of  data  points  on  either  side  of  zero 

NP  - Common  storage  location  for  value  of  LL 

R - Temporary  storage  used  in  calculations 

RJ  - Local  variable  used  in  calculation  of  filter  weights 

RM2  - Storage  location  for  value  of  R-2 

RNR  - Storage  location  for  16  times  deflection  variance 

RRR  - Storage  location  for  2^1  - R + R‘ 

RS  - Storage  location  for  Rc 

Rl  - Storage  location  for  .5[R  -2+2  (1-R+  R“)^]^ 

R1A  - T mporary  storage  of  intermediate  calculation 

RlRJ  - Temporary  storage  used  in  calculation  of  filter  weights 

R2  - Storage  location  for  ,5[-(R-2)  + 2(1-R  + Rs)2]^ 

R2RJ  - Temporary  storage  used  in  calculation  of  filter  weights 
R6  - Storage  location  for  R/6 

-3 

SDD  - Storage  location  for  1 X 10  times  single  pulse  radar  standard 
deviation 

SDDS  - Storage  location  for  square  of  SDD 

SDV  - Storage  location  for  4.84136811  X 10  ^ times  standard  deviation  of 
deflection 
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SDVS  - Storage  location  for  square  of  SDV 

SJ  - Temporary  storage  used  in  calculation  of  filter  weight 

SRA  - Storage  location  for  (1.-  R - R") 

% 

SRB  - Storage  location  for  (1.  + R) 

SRBl  - Temporary  storage  used  in  calculation  of  filter  weights 
SR3  - Storage  location  for  the  constant  1.732050803 
WN  - Array  storage  for  geoid  height  filter  weights 


A 


M 


CMPSGH 

I - Do  loop  variable 

IA  - Variable  used  to  sequence  the  filter  points 
IL  - Number  of  filter  fit  points  plus  one 
J - Variable  counter 
K - Do  loop  variable 

NPT  - Variable  for  total  number  of  data  points  minus  the  total  number  of 
filter  fit  points 

NT  - Variable  for  total  number  of  data  points  minus  % number  of  filter  fit 
points 

NZ  - NPT 

YNN  - Summation  variable  for  geoid  height  or  vertical  deflection  calculations 
YNT  - Array  for  calculated  geoid  heights  or  vertical  deflections 


t 

: 

I 
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1 

. 


V 

' 


i 


STATC0M 


AVE  - Local  variable  used  for  storage  of  average  geoid  height  or  average 
vertical  deflection  residual. 

AXDS1  - Storage  location  for  the  maximum  value  of  geoid  height  or  vertical 
deflection  residual 

AXDS2  - Storage  location  for  the  minimum  value  of  geoid  height  or  vertical 
deflection  residual 

AX1  - Storage  location  for  the  absolute  value  of  AXDS1 
AX2  - Storage  location  for  the  absolute  value  of  AXDS2 
M - Storage  location  of  1,  used  in  PL0T  argument  list 

SS  - Summation  variable  used  to  compute  standard  deviation  of  geoid  height 
or  vertical  deflection  residual 

WN  - Array  used  to  store  geoid  height  or  vertical  deflection  residuals 


DVCMPS 


Definition  of  variables  have  same  meaning  as  used  in  GHCMPTS  except 
that  deflection  of  vertical  calculations  are  performed. 


'I 


PL0T 

DA  - Storage  location  for  geological  feature  being  processed 
DX  - Storage  location  of  the  scaled  minimum  value  for  the  graph  abscissa 

DY  - Storage  location  of  the  scaled  maximum  value  for  the  graph  ordinate 

H - Storage  location  of  the  film  identification 
I - Do  loop  variable 

IP  - Storage  location  for  the  plot  interval 

IT  - Storage  location  of  -1,  used  in  GRF  argument  list 

LA  - Storage  location  of  IP  value,  used  in  GRF  argument  list 

LK  - Storage  location  of  IL  value,  used  in  GRF  argument  list 

Ut  - Storage  location  of  IP  value,  used  in  GRF  argument  list 

LTT  - Local  variable,  storage  for  number  of  points  in  ZN  array 

M - Storage  location  for  0 or  1,  used  in  GRF  argument  list 

ML  - Storage  location  for  the  value  of  the  number  of  points  to  be  plotted, 
used  in  GRF  argument  list 

MP  - Same  use  as  ML  above 

R - Storage  location  for  the  exponent  of  the  scale  factor  which  appears 
on  graph 

TN  - Storage  location  for  Gaussian  constant  which  appears  on  graph 

XHI  - Storage  location  for  maximum  abscissa  object  size 

XL <b  - Storage  location  for  minimum  abscissa  object  size 

XMA  - Storage  location  for  maximum  value  of  graph  abscissa,  used  in  call 

to  GRF  argument  list 

XMIN  - Minimum  value  of  graph  abscissa,  used  in  call  to  GRF  argument  list 
X0R 
X10 
YHI 
YL  <t> 

YMA 
YMX 

YMl 

I 
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Same  as  for  XMIN  above 

Storage  location  for  Hollerith  data  which  goes  on  graph 

Storage  location  of  maximum  graphical  object  size 

Storage  location  of  minimum  graphical  object  size 

Storage  location  at  maximum  ordinate,  used  in  GRF  argument  list 

Storage  location  of  maximum  value  of  residual,  used  in  GRF  argument 
list 

Storage  location  for  the  ordinate  absolute  minimum 


YM2  - Storage  location  for  the  ordinate  absolute  maximum 

YdR  - Ordinate  value  of  graphical  origin  (0);  used  in  GRF  argument  list 

ZW  - Array  containing  the  scaled  residuals  to  be  plotted 
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Appendix  B 

Program  Macro-Flowcharts 


id 


RADAP 


RADA  P 


TAPE  2 

v y 


CALL 

TL  1ST 

CALL 

GHCMPS 

CALL 

CMPSGH 

CALL 

STATCOM 

CALL 

DVCM PS 

CALL 

CMPSGH 

CALL 

STATCOM 

REWIND  2 


SI 


GHCMPTS 


DO  400  I = 1 ,NP 


COMPUTE  SUM  OF  SQUARE 
OF  RESIDUALS  « 
AVERAGE  RESIDUAL 


COMPUTE  STANDARD 
DEVIATION 


STATCOM 


STATCOM 


DVCMPS 


r 

PLOT 


READ  GRAPHICAL 
DATA 


PLOT 


ADD  LABELS  FOR  GRAPH 
WEIGHTING  FUNCTION 
FOR  DEFLECTION  OF 
VERTICAL 


ADD  LABELS  FOR  GRAPH- 
TRUE  GEOID  HEIGHT  AND 

GEOID  HEIGHT  RESIDUALS 


ADD  LABELS  FOR  GRAPH • 

TRUE  DEFLECTION  OF 
VERTICAL  AND  DEFLECTION 

OF  THE  VERTICAL 
RESIDUALS 
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PLOT 


Y/  10. 


o o 


KAO  AP  PROGRAM 


♦G-CR  RAOAP 

PROGRAM  RAiiAF<  INPUT .OUTPUT  ,T APE1,  T 1PE10 , T APE51=OUTPUT , T APE2 > 
COMM ON/ ARRAY/  X<6300),Y<c30Q>,YN<63DO),Dt6100),WN(6300), 

1 YNT  { 4300) 

C ***  CARU  INPUT  *** 

COfoMON/LCIN/  8R0C,  SOO,  SO*,  FOS,  LSX,LFSP 
COMMON/VAR  / NR, DATE,  NG3  , SN  , N 

1,IA,AP,L4. 

COrMON/LOC/L  ,NFS 
CO  MON/NAM:/ H(  2)  ,2  <20  0) 

DIMENSION  *(3) 

PRINT  1000 
PRINT  1010 
PRINT  1020 

***  INPUT  CARD  DATA  *** 

***  I T=2  OR  3 <TYPE  GECICAH.  FEATURE)  *** 

READ  1030,  iRDC  , SOD, NFS  ,LSX,IC 
PRINT  1035,  URuC  ,SDO,NFS  ,LSX,IC 
READ  1040  , SDV  ,FOS,LF3P 
PRINT  1045,  SOW  ,FOS,lFSP 
lFS  =2  * LF5P  ♦ 1 

L =CF  SP 

IGF  =1 

IT  =2 

*EhINU  1 
READ  (1)  NFS 
mENINO  1 
REWIND  2 

IF  (NFS.cU.t#  IT  =3 
IF  ( I T • EG *2  ) GO  TO  200 

READ  (1)  NFS, NR, DATE,  IBTI  ),I=IGF, IT) ,NGS,SN,N 
PRINT  1050 

PRINT  10  55  , NFS, NR, CATE,  <8  < I > , 1= IGF , I T ) ,NGS,SN,N 
IF  (IC.6U.0)  GO  TO  250 
PRINT  1070 
C ♦♦♦  1 APE2- READER  1 *** 

WRITE  (2)  NFS,NR,OATE, <8(1 ), I=IGF , I T » ,NGS , SN , N 
PRINT  1055,  NFS,  NK,  CATE,  (ifi(JI),I  = IGF,IT)  ,NGS,SN,N 
C ***  TAPE2-HEADER  2 *** 

II  = N - 2 * LFSP 

WRITE  (2)  NFS, dRUC.SOO.LSX  ,SOV ,FOS , LFS, 1 1 
PRINT  1075  , NFS»8R0C»SQD, LSI, SDV, FDS, LFS, II 

bO  TO  250 

200  READ  <1)  NFS, NR, DATE, (8(1  l»I=IGF,IT),NGS,SN,N 
PRINT  1050 
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PRINT  lOoO .NFS, NR, CATE, (B ( I > , 1= IGF , I T ) , NGS , SN, N 
IF  (IC.ta.O)  GO  TO  250 
PRINT  1070 

G TAPE2-HEA0ER  1 *** 

HRITE  (2)  NFS, NR  , 0 ATE  , (d(I  ),  1=  IGF  , IT  ),NGS,SN,N 
PRINT  iOcO  , NFS, NR, CATE,  (d C I > , 1= IGF « IT ) ,NGS, SN, N 
C ***  IAPE2-HEA0ER  2 *** 

II  = N - 2 * LFSP 

NRJTt  (2)  NFS, BROC  , S JO  » LS  X ,SDV,FOS,LFS,II 
PRINT  1075  , NFS,BRDC,SOD, LSX, SDV , FDS, LFS ,11 
250  CALL  TLIST  (X  , Y , YN  , 0 ,N) 

CALL  GHCMPTS 

CALL  CMPSGH  ( HN, YNT.NP, YN, 1> 

CALL  STATCOM  <1,  IC) 

CALL  OVCMPo 

CALL  CMPSGH  (WN,YNT,NP,YN,Z) 

CALL  STATCUM  (2,  IC) 

BEHIND  2 

CALL  CKTIQ  (Z,H) 

CALL  EXIT  G (Z) 

10  00  i-  ORMA  T (1 H , * RA  JAR  ALTIMETER  OATA  ANALYSIS  PROGRAM  (RAOAP)  .VERSION 
1-2,11/14/72*,/) 

1010  FORMAT (1H1,41X**RACAR  ALTIMETER  OATA  ANALYSIS  PROGRAM*,//) 

1020  FORMA  T (1H  ,*INPUT  CARO  OATA*,//) 
lOiJ  FORMAT (2t25. 14, 3151 

10  35  FORMA  T ( 1H  ,*BRJC  = *,E25.H,*  ( /KM)  - INVERSE  CORRELATION  OISTANC 
It*,/1H  , 

l*SUO  = *,E25. 14,* (METERS)  - SINGLE  PULSE  RAO AR  SIGMA*, /1H  , 

2*  NFS  = * , 1 5 , * FUNCTION  OF  VERTICAL* , /1H  , 

2*LSX  = * , 15,*  FILTER  FIT  / OATA  SPACING* , /1H  , 

4*  IC  = *,  1 5 , * ( 0-  NC  FILE  , 1=  OUBPUT  FILE)*,/) 

1040  FORMAT (2t2  -.14,110) 

1045  FORMAT ( 1H  ,*SOV  = *,E2I.14,*  (ARCSEC)  - SlfiMA  OEFL ECT ION* , /1H  , 
1*F0S  = *,£25.14,*  (KMETERSI  - OATA  SPACING*, /1H  , 

2*  LFSP  = *,1110,*  1/2  (NUMBER  OF  FILTER  FIT  POINTS)*,/) 

1050  FORMAT ( 1H  ,*TAPE  INPUT*,//* 

1055  FORMA  T (1H  , 215 , 4 A 1 0, 1 5 , E 22 . 1 4, 1 5 , // ) 

10  oJ  FORM  A 1 1 1 rt  , 215, 3A10, 15, E22. 14, 15,//) 

10  70  FORMA  T (1 H ,*TAPE  2 OUTPUT  HEAOER  RECORUS* , // ) 

1075  FORMAT  (1H  , 15 , 2E2 5. 14, 15 , 2E25. 14 , 215 ,// ) 

END 
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......... ...... 
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GHCMPTS 
SUBROUTINE  GHCMPTS 

COMMON/ ARRAY/  X(b300),Y(e300),YN(630Q),Dt6300>,WNC6300), 
1 YNI ( €300) 

***  CAPO  INPUT  *** 

COMMON/ LOIN/  BROC , SOO,  SDV,  FOS,  LSX,LFSP 

9-  9 * . w w 

' " “ . . NGS  p SN  , N 

1, 


TAPE  : 

INPUT  *** 

COMMON/ VAR  / NR, GATE. 

, IL,NP,LL 

COMMON/LOC/L  , NFS 

***  COMPUTE  filter  h 

10 

=L  * U.SX 

LM 

=«N-1)  / 2 

IL 

=LM  - 10  ♦ 1 

SR3 

=1 . 7 32050806 

iOV 

= 4-8A8136611 

sows 

= SOW  * SOW 

SOO 

= l.t-03  * SOI 

SOOS 

= SO u * sou 

dRJCS 

= dROC*  BROC 

oROCU 

= BROC*  dROCS 

RNR 

= 16.*  SOWS 

OR 

= SOOS  * FOS 

R 

= RNR  / OR 

R 

= A6S  (R) 

R 

- ALOG 

R 

= EXP  (R) 

RS 

= R * R 

KM2 

= R - 2. 

SR  A 

= 1 - R ♦ RS 

SRA 

=SURT  (SRA) 

RRR 

= 2 * SRA 

AR 

= 1.  * R 

SRB 

= SORT  (AR) 

R1A 

= RRR  ♦ RM2 

R1 

= SORT  (R1A) 

R1A 

= RRR  - RM2 

R2 

= SORT  (R1A) 

Rfe 

= R / 6. 

LL 

= 2 ♦ LFSP  ♦ 

NP 

=LL 

J 

=-  lFSP  - 1 

GO  2 

00  i = 1,  LL 

J = J 

♦ 1 

FDS1 

=J  * FOS 

S J 

= AuS  ( FOS1) 

RJ 

= BROC  * SJ 

R1RJ 

= R1  * RJ 

R2RJ 

= C -R2 ) * RJ 

oRBl 

- - 6RB  • RJ 

E-06  * SDV 


BROCfl 


(R)  / 3. 


.5 


.5 
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s. 


"Vi  • : 


F = R6  * l (Rl*  (SIN  (R1RJ)  ♦ SR3  MCOS  IR1RJ))»  ♦ 

1 R2*  ( SR3*SIMR1RJ> - <COS  (R1RJ)»>>* 

2 SR A ♦ EXP  (SRB1 ) / SRB> 

WN  (I)  = BROC  ♦ F 

200  CONTINUE 
PRINT  400 

CAUL  POUMP  (WNT1I » WN  (LIL )* 1 1 

CALL  PLOT  ( WN  , X ,LL,  1»  1 * NFS,IL) 

*00  FORMAT  (1H  ,*GEOIO  HEIGHT  WEIGHT* »//l 
RETURN 
ENO 


EXP (R2RJ ) 
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* KECK  CMPgGH 

5U6R0U1INE  CJfPSGH  (WN,YNf!,NP,YN,IC) 

G ”*  gamD  input  *♦* 

COFMON/LCIN/  broc,  sjo,  sdv,  fos,  lsx.lfsp 
COPMON/VAR  / NR, DATE,  NGS  , SN  , N 

1, IL,NZ,LL 
COPMON/LOC/L  ,NFS 

DIMENSION  rNirC6300),WN<Ni»)  , YN(6300) 

IL  = (ifcP  - 1>  / 2.  + 1 

NPT  =N  - 2.  * (IL  - 1) 

NT  = IL  + NPT  - 1 

NZ  =NP T 

J =0 

00  210  1 = IL  » N T 

J = J +1 

YNN  = 0. 

DO  200  K = 1 » LL 

IA  = J ♦ (K  - 1)  * L3X 

YN#t  = YNN  ♦ WN  ( K)  * YN(IA) 

200  CONTINUE 

YNT  (J  ) = FOi  * YNN 

210  CONTINUE 

IF  (IC.EQ.2)  GO  TO  400 
PRINT  300 

300  FORMAT  (1H  ,*GEOIO  HEIGHT  fSTIMATES*,//) 

GO  TO  500 
*♦00  PRINT  350 

350  FORMAT  <1H  , ’VERTICAL  DEFLECTION  ESTIMATES’,//) 
500  CALL  PUUMP (YNT (1) , YNT (NP» , 1) 

RETURN 

END 
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DECK  STATCOM 

SUBROUTINE  STATCOM  (IC,IIC) 

COMMON/ ARRAY/  X(6300)  ,Y  (6300  ) ,YN(t63  0 0)  ,0  ( 6300  ) , UK  ( 630  0 ) , 
1YNT( E3U0) 

CC*NON/LCIN/  BRDC,  SOD,  SOW,  FDS,  IISX.LFSP 
COMMON/ VAR  ✓ NR, DATE,  NGS  , SN  , N 

1,  IlL,  NP,LL 
COMMON/LOC/L  ,NFS 

***  COMPUTE  RESIDUALS  OF  EST  ANO  TRUE  VALUES  *** 

IF  (IC.EU.2)  GO  TO  150 

UO  100  I = 1,  NP 

WN  (I)  = YMT  ( I > - Y (I  ♦ IL  > 

100  CONTINUE 

IF  (IIC.EQ.U)  GO  TC  300 
WRITE  (2)  NP, ( WN  C I ) , I = 1,NP> 

GO  TO  300 

150  UO  200  I = 1,  NP 

WN  (I ) = YNT  « X > - J Cl  ♦ IL) 

200  CONTINUE 

IF  (IIC.EQ.U)  GO  TO  300 
WRITE  (2)  NP,  ( HN  ( I ) , I = 1,NP) 

300  AXOS1  = AMAXE  (WN  , NP) 

AXUS2  = AMINE  IWN  , NP) 

AX1  = AtiS  ( AXOS1) 

AX2  = AiS  ( AX0S2) 

A XUS  1 = AMAX1  ( A XI  , «X2) 

COMPUTE  AVE  VALUE  OF  RESIDUALS  *** 

A V£  = 0. 

UO  350  I = 1,  NP 
AVE  = AVt  WN  (I ) 

350  CONTINUE 

AVE  = AVE  / NP 

COMPUTE  STANOAFO  DEVIATION  OF  RESIDUAL 
SIGMA  = 0. 

UO  400  I = 1,  KP 
SS  = WN  (I)  - A Vi 

SS  = S3  * SS 

SIGMA  = SIGMA  ♦ SS 
400  CONTINUE 

SIGMA  = SIGMA  / NF 
SIGMA  = jukT  (SIGMA) 

IF  (IC.tQ.2)  GO  TO  450 
PRINK  1000,  AXOS1 
PRINT  1005„  AVE 
PRINT  1010,  SIGMA 
M =1 

CALL  PLOT  I Y , X , N,  3,  2 , NFS , M ) 

RETURN 

450  PRINT  1015  , AXOS1 
PRINT  1020  , A VE 


FRlNl  1025  ,SIGMA 
M =1 


CALL  PLOT 

1 0 

* x » 

N,  «♦, 

2 , NFS , M ) 

10  00 

FORMAT (1H 

, *MA 

XI  MUM 

GEOIC 

HEIGHT 

RESIDUAL  = *, 

Ffc.2,* 

METERS* 

»/) 

10  0 3 

FORMA  T ( 1H 

,*  AV 

tRAGE 

uEOIB 

FEUGHT 

RESIOUAL  = ♦, 

Ftt.  2,* 

METERS* 

,/) 

1010 

FORMA  T ( 1H 

,*SI 

GMA  CF 

GEOID 

HEIGHT 

ESTIMATE  = ♦ 

,F6.2  ,* 

METERS* 

,/) 

1015  FORMAT  (1H  ,*MAXIMUM  VERTICAL  DEFLECTION  RESIOUAL  = *,F6.2,*  ARCSEC 

1*,/) 

1020  FORMA  T <1 h ,*AVERAGE  VERTICAL  DEFLECTION  RESIOUAL  = * 
lu201,F€.2,*  ARc;>EC*,/> 

1025  FORMAT ( 1 H ,*SIGMA  OF  VERTICAL  DEFLECTION  ESTIMATE  = * 

10251,F6.2,*  ARCSEC*,/) 

RETURN 

cND 
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♦GcCK  JVCMPS 

SUBROUTINE  DVCMPS 

LOrtMON/A*RA V/  X(6300),Y(fc300>,YN(6300>,Ot63aO>,WN(630G), 

1YNT ( 43U0) 

C ***  CARO  INPUT  *♦* 

C 0 R M 0 N/  L C I h / CROC,  SOLI,  SUV,  FOS,  ILSX,LFSP 
C ***  TAPE  INPUT  *** 

COMMGN/VAR  / NR , 1) A TE , NGS  , SN  , N 

1 , I L , N P , LL 
COPNON/LOC/L  , NFS 

C ♦**  CCMPUTc  FILTER  WEIGHTS  *** 

j.0  =L  * dSX 

LM  = (ih  - 1)  / 2 

IL  =LM  - 10  ♦ 1 

SR3  =1. 732050  80  E 

SGVS  = SGV  * SUV 

SOUS  = SJO  * soo 

URUCS  = BRUC*  8RGC 

BRQCU  = BROC*  BROCS 

RNR  = 16.*  SUVS 

uR  = SOUS  * FOS  * BRUCG 

R = RNR  / OR 

R = A2S  (R) 

R = ( ALGG(R) )/ 3. 

R = EXP  ( R I 

RS  = R * R 

RM2  = R - 2. 

BRA  = 1 - R * RS 

RRR  = 2 * SURT  (SR A) 

AR  =1.  ♦ R 

SRti  = SURT  (AR) 

RlA  = RRR  * RM2 

R1  = SURT  (RlA)  * .5 

RlA  = RRR  - RM2 

R2  = SURT  (RlA)  * .5 

R6  = R / 6 . ' 

LL  = 2 * LFSP  ♦ 1 

NP  =LL 

J =-  lFSP  - 1 

C CALL  PGUMP(R*R,RNR,1,R,R,1,R1,R1,1,R2,R2,1)) 

UO  2 00  1 = 1,  LL 

J = J ♦ 1 
FOS1  =J  * FOS 

BJ  = AoS  ( FOS1) 

RJ  = BROC  * SJ 

R1RJ  = R1  * RJ 

k2RJ  = ( -R2 ) * RJ 

SR3E  =-SRB  ♦ RJ 

F = R6  * ((  COS(RIRJ)  - SR3  * SI*  (R1RJ) ) * EXP  (R2RJ) 

1 - EXP  (SR3E ) ) 
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IF  CFDii. EC. 0.)  GO  TO  190 
WN  Cl)  = BKOCS  * F*FOSl  / SJ 
GO  TO  200 

190  WN  (I)  = 3RUCS  * F 
200  CONTINUE 
PRIN1  400 

«*00  FORMAT  (1H  , ♦VERTICAL.  REFLECTION  W £1  (6HT1S*  , // ) 
CALL  FuUMPC  WN( 1» ,WNCNP),1,IL,IL,2) 

CALL  FLUT  ( WN  , X ,LL,  2,  1 , NFS  * IL I 

00  300  I = 1 » LL 

WNCI)  = 206.2648062  * HI«(I) 

300  CONTINUE 
RETURN 
ENO 


O o 


I, 


♦O^CK  PLOT 

SU0RGUTINE  PLOT  CZN,X  ,NL  , LT ,L T 1 , NFS, LK) 
COMMON/LOIN/  6R0C,  SOJ,  iUL,  FQS,  I1SX,LFSP 
COMMON/VAR  / NR.OATE,  NGS  , SN  , N 

1,IL,NP,Ll 
COMMON/NAMc./  H ,Z 

COMMON/ ARRAY/  OUM ( 25 2 00 ) »WN C 630 U) ,OUM1(6300> 
UIMENS10NZN  ( -6.3 □ 0 ) , X(NL)  ,ZH  16300  ) 

UIMENSION  M2  > , Z 1 200 ) 

UIMLASiON  G A ( 3 ) , TM3  ) 

U AT  A XLO, YLO.XHI, YHI  / . 0 3 , . 0 6 , 1 . 14 , .94/ 

(DATA  X10/3PX10/ 

LT  = G*APH  NO. 

***  L T 1 = l OR  2 PLOTS  / GRAPH 
CALL  MGbcSC  (Z • Q > 

IF  (LT.GT.U)  GO  TO  100 


READ 

1000  , H 

PRINT 

1005  , H 

CALL 

CkTIO  (Z,H) 

READ 

1010,  OA , T N , IP 

PRINT 

1015  ,DA,TN,IP 

loo  ltt 

= nl 

YMI 

= AMINE  1ZN,  LTT  ) 

YMX 

= AMAXE  ( ZN , LTT  > 

YMI 

= A&S  1 YMI ) 

YM2 

= A3S  ( YM  X 1 

YMI 

= AMAX1  ( YM 1 , Y M2 ) 

IFJLT1.EQ.1)  GO  TO  150 


YMI 

= 

AMINE 

(HN,NP> 

YMX 

= 

Ari  AXE 

( WN ,NP ) 

YMI 

= 

A8S  (YMI) 

YM2 

= 

AOS  (YMX) 

YMI 

= 

ArtAXl 

(YMI, YM2) 

YMA 

= 

AMAX1 

(YMI, YMI) 

XMIN 

X ( LK) 

XMA 

= 

-XMIN 

CALL  SCALE  (YMA  ,KE,YMI  I 
YM A = -YMl 

XOR  = X N 1 N 

YOiR  = 0. 

OX  = 2 * XliA  / 10. 

L» Y =Y  M A/5  . 

R =-Ks- 

IFlLll.cQ.l)  GO  TO  205 
00  2U0  I = l,NP 

HN  (I)  = WN(I)  * 10.**  R 
2 0 J CONTINUE 
205  00  210  1 = 1 , NL 

ZH  (I ) = Z N ( I > * 10.**R 


C-ll 


210  COM  I NUt 


IT 
t.  A 
lO 
M 


= -i 

= IP 
= IP 
= 0 


3 0 13  CALL  SUBJEG  ( Z , 0 . , C.  , 4095  . ,3  0 71  . > 
CALL  OBJCTG  HZ,XLU,  YLO.XHl,  Y HI  ) 
CALL  Stl 3MG  (Z,  45,  .75) 

CALL  ScTSMiI  ( Z , 1 1 3 ,XLO) 

CALL  StTSMG  «Z,114,Yl0) 

CALL  StTSMG  CZ,115,XHI) 

CALL  SlTSMl  ( Z , life  , YH I ) 

CALL  NuMdRC  (Z , 02 4 8.  ,2 5 3 1. >3  ,KE) 
CALL  LlGNDG  < Z , 0 1 8 6. , 2h79 . ,3  ,3HX10) 
bO  TO  (350,450,550,65  0)  ,L  1 


350  CALL 

LlGNOG 

(Z  ,1561. ,273b. 

, 35, 

35HHEIGHT ING  FUNCTION  FOR 

1GHT) 

CALL 

LLGftOb 

<Z , 1 56 1.  ,2683. 

, 9, 

9HFUN  NO.  ) 

CALL 

NUMdRG 

(Z,  1840. ,2683. 

, 5. 

NR  ) 

CALL 

LLGNUG 

(Z,  2150. ,2683. 

, 6, 

6H0ATE  ) 

CALL 

NUMbRb 

(Z,2460.,2683. 

,-io. 

DATE  ) 

CALL 

L c.  GU  0 G 

IZ,  1561. , 2632. 

, 15, 

15  HO AT  A/N  OISE  TYPE) 

CALL 

Lt&NOG 

(Z  ,1933. ,26 32. 

, 30, 

OA) 

CALL 

LlGNOG 

(Z , 2 52  6.,  2632. 

, 3, 

3 H / ) 

call 

LlGNOG 

(<Z  , 1 9 3 3 • ,2581. 

, 30, 

TN) 

CALL 

TITLEo 

<Z,0,XC,14,14HFILTER 

HEIGHTS, 0, TC) 

call 

TiUEiUZ,23,23HAL0NC 

1RACK 

UISTANCE-KM,  0,YC,0  ,TC) 

GO  TO 

7 00 

+50  CALL 

LlGNOG 

(Z, 1561. ,2734. 

» 45, 

45HWEIGHTING  FUNCTION  FOR 

ION  OF 

VERTICAL) 

CALL 

Lt  GNU  C 

(Z,1561. ,2683. 

, 9, 

9HRUN  NO.  ) 

CALL 

NOMtiRL 

(Z  ,1840. , 268i. 

. 5, 

NR  ) 

CALL 

LlGNOG 

(Z  ,2150. , 2683. 

, 6, 

&HDATE  > 

call 

NUMdRI 

<(Z  ,2460.  ,2683. 

,-io, 

OATE  ) 

lALL 

Li_GN  Ob 

«Z  , 1 5 6 1 . ,2632. 

, 15, 

15HU AT  A/NOISE  TYPE) 

call 

Lt&NUG 

(Z, 1933. ,2632. 

» 30, 

□ A) 

lALL 

LcGNuG 

(Z,2526. *26 32. 

, 3 , 

3 H / ) 

CALL 

LlGNOG 

HZ,  1933. ,2581. 

, 30, 

TN) 

CALL 

UTlEG(Z,U,XC,14,14HiFILTER 

wei  Gifts  , o , no 

CALL 

TITLEO 

(Z, 23, 23HAL0NG 

TRACK 

U I ST  ANCE -KM , U,YC,0  ,TC> 

bO  10 

700 

550  bALL 

LLONi)o(^yl361«)27i4«  9 

21,21HTRUE  GEO  1 0 HEIGHT  ANO) 

CALL 

L_GNUG(Zf , 

12 , 22HGE0 10  HEIGHT  RESIDUALS) 

CALL 

LiGhUo 

CZ,15bl. ,2632. 

, 9, 

9HRUN  NO.  ) 

CALL 

NUMUR  i 

IZ,  1340. ,2632. 

, 5, 

NR  ) 

CAIL 

L-GNOG 

IZ  ,2150. ,2632. 

, 6 , 

6FDATE  ) 

CALL 

NuMdKo 

(Z  ,2460. , 2632. 

,-10, 

U ATE  ) 

CAU.L 

LlGNOG 

(Z,  1561.  ,25811. 

, 15, 

15H0ATA/N0ISE  TYPE) 

CALL 

L-GNOb 

«Z, 1933. ,2581. 

, 30  , 

OA) 

b AlL 

LEGNGb 

IZ  ,2526. ,2581. 

, 3 , 

3H  / ) 

o r.  o o 


CALL 

LlGNuG 

(Z,  1933. ,2531. 

, 30,  TN) 

CALL 

TITL£6(Z,U,XC  ,19,19HGE0I0  HEIGHT-METERS 

, 0 , TC ) 

CALL 

TiTlEMZ ,23, 2 3H A LONG 

IRACK  OISTANCE-KM, 

0 ,YC,  0 

,TC) 

GO  TO 

700 

o50  CALL 

LcGNL)  G 

(Z  , 1561. , 2 73h. 

,3 1 , 31H  TRUE  DEFLECTION  OF 

VERTICAL 

CA  X 

LEGNUo 

(Z , 1561. „2683. 

,3fc,36HDEFLECTI0F 

OF  THE 

VERTICAL 

1 ALS ) 

CALL 

LlGNDG 

CZ, 1561. ,2632. 

, 9,  9PRUN  NO. 

) 

CALL 

NUMiIRl 

tZ , 1 ft A 0 . , 2 6 32  . 

, 5,  NR 

) 

CALL 

L-GNOG 

OZ, 2150. ,2632. 

, 6,  6HDATE  ) 

CALL 

NCHoRvb 

U ,2460. , 2632. 

,-10,  DATE  ) 

CALL 

LEGNOG 

«Z, 1561. ,2581. 

, 15, 15H0AT A/NOISE 

TYPE) 

CALL 

legnog 

«Z, 1933. ,2581. 

, 30,  OA) 

CALL 

LlGNuG 

OZ, 2526. ,2581. 

, 3,  3H  / ) 

CALL 

L-GNOG 

(Z,1933.,253i. 

, 30,  TN) 

CALL 

TlTLcG(Z,Q,XC,29,29H0EFLECTI0F  OF  VER T I CAL- ARCSEC ,0 ,TC) 

CALL 

TITLEGTZ, 23, 2 3H ALONG 

TRACK  0 1 ST ANCE-Kt* , 

0 , YC , 0 

, TC) 

700  CALL 

GETSMG 

0Z,40,Z4Q> 

CALL 

Gc.T  SMG 

ttZ , A 1 , Z41 ) 

CALL 

CoJC  T G (Z , XLO  ♦Z.’ZAC- 

2 , * ZmO  , 

1YL0*2.*ZA1 ,XHI-2.*Z40 , YHI-10.*Z41> 

CALL  SclSMG  <Z,55,0  .) 

CALL  POUHP  CXMIN,XMIN,1, Yfl  ,Y M I , 1 , X X A ,XMA  ,1,YMA ,YMA  , 1 , XOR ,X OR , 1 , YOR 
l,Y0R,l,0X,fcX,l,0Y,0Y,l,NL,NL,2,X<LK>  , X <NL  > , 1 , ZW  ( 1)  , ZM(NL),  1) 

ML  =NL/IP 

CALL  GkFUXMIK,  YMI  , XM  A,  Y If  A , XOR  , YOR  ,OX , D Y , 

10. 2.0,  ML,  X ( LK ) , ZH  (1),  LA  , LO  , 1H*  , 

21 T , M i ll  • i 0 » 0.  j X l 1 • 1 • f A , Of  Z) 

IFTLTl.tO.l)  GO  TO  750 

LK  =IL 

CALL  S^TiMG(Z, 55,0.) 

CALL  FoUMP<lL,IL,2,KE,KE,2,YMA, YMA.l) 

CALL  PJU.1P  <NF,NP,2  ,HN  ( 1 > »Wit  (NP  > ,1  ) 

MP  =NP//1P 

CALL  GKFIXMN,  YMI  , XM  A , YPA  , XOR  , YOR  ,OX  , OY  , 

10.2.1, nP,  X ( L K ) , MN  (1),  LA  , LO  , 1H  + , 

2IT  , M , i)«  , 0 , 0.  , X L 1 , 1.  , h , 0,  Z) 

750  CONTINUE 

CALL  PAGcGlZ, 0,0,1) 

1000  fORMAT  (2  A 10  ) 

10  05  FORMAT (1H  ,*FILM  I CENT  IF IC AT  ION  = *,2A10) 

1010  FORMAT  C3A1U, 3A10, 15) 

1015  FORMA  M1H  ,*OA  = *,3A10,*TN  = *,3A10,*IP  = *,15) 

RETURN 

END 


Si 


C-13 


*OECK  SCALE 

SUBROUTINE  SCALE  <YMA,I  ,YH  ) 

C ***  IMA  GT  1 *** 

I =0 

IF  (YMA.LE.l.)  GO  TO  500 

Y = YM  A 

100  I =1*1 

Y = Y ✓ 10 

IF  (T.GT.l.)  GO  TO  100 

C ***  YMA  IN  FRACTION  FORM  .IMA  X 10  >*  I) 

IF  (Y.GT..5)  GO  TO  150 
YMI  = -.5 

RETURN 

150  YMI  =-l. 

RETURN 

500  IF  (YMA.EQ.l.)RETURN 

Y = YMA 

550  I = I *1 

r = la  * y 

IF  CY.LT.  1.)  GO  TO  550 
I =1-1 

Y = Y / 10. 

I = -I 

IF  (Y.GT..5*  GO  TO  600 

YMI  = -.5 

RETURN 

600  YMI  =-i. 

RETURN 

ENO 


C-14 


Appendix  D 

Sample  of  Program  Print  Out 
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HI-D  BLAKE  ESCARPMENT  401CA  TN  = GAUSSIAN  SIGMA=  0.6  M 


c 

w 

O 

3 

o 

c/i 

r 

5 

5 

1 

bJ 

o 

Ufir-KTINS  FONCTION  FOB  BFOIO  HEIGHT 
F'.'N  NO.  «M  Oftlf  “>  J»N  '3 

COT  ft  / NOISE  TYPE  NOMNftl  fiEOIO  / 

r,l"r.  5 1 AN  5 I G«A  1 . 3 "I 


ANOF -«« 


AND 
'UAL  S 

DATE  T Al 
NOMINAL  GEOID  2 0(/ 
GAUbblAN  blGMA  = 


r,H3  J 3W-  J My  J 3"  0/033 
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